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Abstract. We present results on the QCD equation of state, obtained with two different improved dynam- 
ical staggered fermion actions and almost physical quark masses. Lattice cut-off effect are discussed in 
detail as results for three different lattice spacings are available now, i.e. results have been obtained on 
lattices with temporal extent of N T = 4, 6 and 8. Furthermore we discuss the Taylor expansion approach 
to non-zero baryon chemical potential and present the isentropic equation of state on lines of constant 
entropy per baryon number. 
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1 Introduction 

A detailed and comprehensive understanding of the ther- 
modynamics of quarks and gluons, e.g. of the equation 
of state is most desirable and of particular importance 
for the phenomenology of relativistic heavy ion collisions. 
Lattice regularized QCD simulations at non-zero temper- 
atures have been shown to be a very successful tool in an- 
alyzing the non-perturbative features of the quark-gluon 
plasma. Driven by both, the exponential growth of the 
computational power of recent super-computer as well as 
by drastic algorithmic improvements one is now able to 
simulate dynamical quarks and gluons on fine lattices with 
almost physical masses. 

In this article we present results on bulk thermody- 
namic quantities and the equation of state on lattices with 
temporal extent of N T = 4, 6 PQ and preliminary results 
on N T — 8, obtained by the HotQCD collaboration[21[3]. 
The article is organized as follows, in Sec. [2] we discuss 
details of our lattice actions and corresponding finite cut- 
off corrections, in Sec [3] we present our choice of lattice 
parameter, in Sec. [3] we present preliminary results on the 
equation of state and in Sec. [5] we introduce the Taylor 
expansion method and calculate leading and next to lead- 
ing order corrections of bulk thermodynamic quantities to 
a non-zero baryon chemical potential. Finally, we present 
our preliminary results on the isentropic equation of state 
in Sec. [5J and conclude in Sec. [7] 



2 Choice of action and cut-off effects 

In order to control lattice cut-off effects it is of particular 
importance for finite temperature calculations to improve 
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the lattice action beyond the naive discretization scheme. 
As most bulk thermodynamic observables, like the energy 
density and pressure, are dimension four operators, the 
numerical signal for these observables drops like the lat- 
tice spacing to the fourth power. One is forced to perform 
lattice calculations on rather coarse lattices, where cut-off 
effects are still sizable. Using an improvement scheme as 
introduced by Symanzik [4] for the gauge part of the action 
and similarly for the fermionic part in the case of staggered 
fermions 0|BJ or Wilson type fermions [3IHEI , cut-off ef- 
fects can be drastically reduced. By adding irrelevant op- 
erators to the action which will vanish in the continuum 
limit, cut-off effect to arbitrary order in 0(a n ) can be 
eliminated already for finite lattice spacing a. The above 
mentioned improvement schemes for staggered fermions 
are tree-level improvements 0(g ) but can in principle be 
generalized to eliminate cut-off effect also in leading order 
of the gauge coupling 0(g 2 ). 

The basis for the two actions we will consider here 
is given by two versions of O(a 2 g ) improved staggered 
fermions which contain quark and anti-quark fields sep- 
arated by up to three links. By adding a straight 3-link 
term to the standard 1-link term, and choosing the coef- 
ficient such that the leading order cut-off effects cancel, 
one arrives at the Naik action [5J. Using bended 3-link 
terms instead (knight moves), one derives the p4 action 
[BJ. The latter action has a dispersion relation which is 
0(p A ) improved. 

In addition to the tree-level improvements one intro- 
duces a smeared I-link term which reduces the flavor sym- 
metry breaking of the staggered fermions. In case of the 
p4-action, the 1-link term is smeared by the sum of all 
corresponding 3-link staples. In case of the Naik action, 
staples up to length of 7 links are used as well as an addi- 
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Fig. 1. Pressure for three different types of free lattice fermions 
as function of N T . Results are given in units of the contin- 
uum ideal Fermi gas value and are calculated for two different 
masses and two different values of chemical potentials. 



tional tadpole improvement (1-loop level). The gauge part 
of both actions is Symanzik improved. Using all this im- 
provements together these actions are called p4fat3 and 
asqtad, respectively. 

In the free gas limit (T — > oo), where the smearing 
becomes irrelevant, the cut-off effects of the p4 and Naik 
actions for bulk thermodynamic quantities have been an- 
alyzed in detail [TO] . Within a large N T expansion, where 
N T is the number of lattice sites in temporal direction one 
can quantify how the pressure approaches its continuum 
(Stcfan-Boltzmann) value psb- In comparison we find for 
the standard (non- improved), p4 and the Naik actions 
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As one can see, the leading order corrections of the im- 
proved p4 and Naik actions are identical. However, the 
sub-leading coefficient which is of order (ir/N T ) 6 is much 
smaller for p4 than for the Naik action. In Fig [T] we show 
the exact evaluation of the pressure as function of N T , 
for two different masses and two different chemical poten- 
tials. It is evident, that the dependence of cut-off effects 
on quark mass and chemical potential /i is very small and 
similar in all three discretization schemes. In fact, one can 
show that the dependence of finite cut-off corrections on 
ji is given by Bernoulli polynomials and is independent of 
the discretization schemes |10) . 



3 Lattice parameter and scale setting 

We perform calculations on lattices of extent 16 3 x 4, 
24 3 x 6, using the p4fat3 action Q] and 32 3 x 8 with both 
p4fat3 and asqtad actions. The latter calculations are still 
preliminary and are currently performed by the HotQCD 
collaboration [2]. For the generation of gauge configura- 
tions we use the exact RHMC algorithm [TT]. For each 
finite temperature calculation we perform a correspond- 
ing zero temperature calculation on a lattice of at least 
the size N 4 , where N„ is the spatial extent of the finite 
temperature calculation. 

The simulations are done on a line of constant physics 
(LCP), i.e. the quark masses are kept constant in phys- 
ical units. In practice, this has been obtained by tuning 
the bare quark masses such that the meson masses of e.g. 
pion, kaon and pseudo-scalar strange meson ss stay con- 
stant in the QCD vacuum as we change the value of the 
coupling. The strange quark mass was always fixed to its 
physical value, by fixing kaon and ss to their correspond- 
ing physical values [TJ. We find that the LCP can, to a 
good approximation, be parameterized by a constant ra- 
tio of the bare quark masses. Most calculation are done 
on a LCP which corresponds to a pion mass of about 220 
MeV. We do, however, also show preliminary results with 
a physical pion mass, i.e. m w « 150 MeV. 

To set set the temperature scale in physical units, we 
determine two distance scales, r$ and ri, from the zero 
temperature static quark potential 



dr 



1.65, 



2 dV qq {r) 
dr 



1.0 (2) 



The ratio of both scales is only slightly quark mass de- 
pendent. It has been determined in both discretization 
schemes consistently, ro/Vi = 1.4636(60) (p4fat3 [TJ) and 
1.474(7)(18) (asqtad [12]). The distance scales r and n 
have been related to properties of the chamonium spec- 
trum which allows to determine them in physical units. 
We use here ro = 0.469(7) fm as determined in Ref. [T3] , 
More details on the scale setting procedure, as well as the 
parameterization of the LCP are given in Ref. [TJ. 



4 The equation of state 

Along the line of constant physics, at sufficiently large vol- 
ume and at zero chemical potential, the temperature is the 
only intensive parameter that controls the thermodynam- 
ics. Consequently there exists only one independent bulk 
thermodynamic observable that needs to be calculated. 
All other quantities are than obtained by using standard 
thermodynamic relations. On the lattice, it is convenient 
to first calculate the trace anomaly in units of the fourth 
power of the temperature, O^/T 4 . It is easily obtained 
as a derivative of the pressure p/T 4 , with respect to the 
temperature, 
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Fig. 2. The trace anomaly, (e — 3p)/T 4 on N T = 6 and 8 lattices. On the left panel results are obtained with the p4fat3 and 
asqtad actions, where N T — 6 results are from [T] and [TJ] respectively, N T — 8 results are preliminary hotQCD results [2]. The 
right panel shows the low temperature part of the trace anomaly in detail, here we plot only results obtained by the p4fat3 
action. Also shown on the right panel are quadratic fits to the data, as well as the trace anomaly obtained in the framework of 
the Resonance gas model. Light quark masses have been constrained to be one tenth of the strange quark mass (0.1m s ), on the 
right panel we show, however, also preliminary results from simulations with physical quark masses (0.05m s ). 



As the pressure is directly given by the partition function, 
p/T = V~ 4 \nZ, the calculation of the trace anomaly re- 
quires only the evaluation of rather simple expectation 
values. According to Eq. [3] one then obtains the pressure 

by 

^-^ = £dT'JL -(T') ■ (4) 

Here Tq is an arbitrary temperature value which usually 
is chosen in the low temperature regime where the pres- 
sure and other thermodynamic quantities are suppressed 
exponentially by Boltzmann factors corresponding to the 
lightest hadronic states; e.g. the pions. The energy den- 
sity is then obtained by combining results for p/T 4 and 
(e — 3p)/T 4 , respectively. 

In Fig. [2] (left) we show results for /T 4 obtained 
with the asqtad and p4fat3 actions, respectively. The new 
N T = 8 results [2] are compared to N T = 6 results taken 
from [TlfTl] . We note that the asqtad and p4fat3 formula- 
tions give results which are in good agreement with each 
other. In fact, in quite a large temperature regime the 
agreement for given lattice extent N T seems to be better 
than one could expect in view of the overall cut-off depen- 
dence that is visible when comparing results for N T — 6 
and N T = 8 more closely. They lead to a reduction of the 
peak height in G^/T 4 , which is located at T ~ 200 MeV 
and to a shift of the rapidly rising part of G^ /T 4 in the 
transition region to smaller values of the temperature. 

In Fig. [5] (right) we show quadratic fits to the data 
obtained by the p4fat3 action, to highlight the cut-off ef- 
fects. It is evident, that the N T — 8 data is shifted relative 
to the N T — 6 data by about 8 MeV at low temperatures, 



T ~ 160 MeV. This shift decreases to about 5 MeV at 
temperatures T ~ 190 MeV. 

Light quark masses have been constrained to be one 
tenth of the strange quark mass (0.1m s ), on the right 
panel we show, however, also preliminary results from sim- 
ulations with physical quark masses (0.05m s ). This again 
leads to a further shift in 6> w of approximately 5 MeV at 
T ~ 190 MeV towards lower temperatures. 

We also compare the results for (e — 3p) /T 4 to results 
obtained from the hadron resonance gas model. Details on 
the resonance gas curve in Fig. [2] (right) will be given in 
[2]. The slope of (e — 3p)/T 4 obtained by the resonance 
gas model seems to be much smaller than the slope ob- 
tained by the quadratic fits to the data. Whether this 
points at deviations of the equation of state at lower tem- 
peratures from resonance gas behavior or is due to larger 
cut-off effects in the low temperature regime requires fur- 
ther studies. We note that the lattice spacing becomes 
larger at lower temperatures and violations of flavor sym- 
metry, which are inherent to the staggered fermion for- 
mulations at finite lattice spacing, thus may become more 
important. 

The cut-off dependence observed in <9 M/i /T 4 carries 
over to the calculation of pressure and energy density; the 
former is obtained by integrating over (9 W /T 5 and the 
energy density is then obtained by combining results for 
p/T 4 and (e — 3p)/T 4 . This is apparent in Fig. [3] where 
we show the ratio p/e obtained with the p4fat3 action 
on N T = 6 [1J and JV T = 8 [| lattices. Cut-off effects 
are still visible in the vicinity of the 'softest point' of the 
EoS, which is related to the peak position of (e — 3p)/T 4 . 
We find that in the entire range of energy densities rele- 
vant for the expansion of dense matter created at RHIC, 
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Fig. 3. The ratio of pressure and energy density as well as 
the velocity of sound obtained in calculations with the p4fat3 
action on N T — 6 [T] and N T = 8 [5] lattices. 



e < 10 GeV/fm 3 , the ratio p/e deviates significantly from 
the conformal, ideal gas value p/e — 1/3. 

This also is reflected in the behavior of the velocity of 
sound, c 2 — dp/de, which is shown in Fig. [3] by dashed 
lines. It starts deviating significantly from the ideal gas 
value below e ~ 10 GeV/fm 3 and reaches a value of 
about 0.1 in the transition region at energy densities e ~ 
1 GeV/fm 3 . Below the transition it slightly rises again, but 
note, that for very small temperatures c 2 , as well as p/e 
are sensitive the integration constant po(T ) (see Eq. [3|). 
At present we have set p (T ) = 0, for T = 100 MeV. 



5 Non-zero chemical potential 

At non-zero chemical potential, lattice QCD is harmed 
by the "sign-problem" , which makes direct lattice calcu- 
lations with standard Monte Carlo techniques at non-zero 
density practically impossible. However, for small values 
of the chemical potential, some methods have been suc- 
cessfully used to extract information on the dependence 
of thermodynamic quantitites on the chemical potential. 
For an overview see, e.g. [T5] . 

We closely follow here the approach and notation used 
in Ref. [16]. We start with a Taylor expansion for the 
pressure in terms of the quark chemical potentials 



P_ 



E 



(5) 



The expansion coefficients c™'?£(T) are computed on the 
lattice at zero chemical potential, using stochastic estima- 
tors. Some details are given in [T7] , 

In Fig. |4] we show results on the diagonal expansion 
coefficients with respect to the up-quark chemical poten- 
tial up to the six order (c^'q'q with n — 2,4,6), obtained 
with the p4fat3 action. Here the full symbols are from 
N T — 4 lattices, while the open symbols denote results 



from N T — 6 lattices. We find that cut-off effects are 
small and of similar magnitude as those found for the trace 
anomaly 6> MM . This was already anticipated by the analy- 
sis of the cut-off corrections in the free gas limit. Similar 
results for the asqtad action have been obtained in [TS] , 

We also compare our preliminary results for 2+1-flavor 
QCD and a pion mass of « 220 MeV, with previously 
obtained result of 2-flavor QCD and w 770 MeV [TB] 
(also p4fat3 and N T = 4). It is apparent from Fig. U]that 
the critical temperature for these two particular sets of 
lattice parameter differ substantially and in fact decreases 
from about 225 MeV for the heavier mass calculations to 
about 200 MeV for the lighter mass calculations. Note, 
that those T c values are the N T = 4 values, which of course 
are still influenced by the finite lattice spacing. Further- 
more, we find from Fig. 3] that the quark number fluctua- 
tions of second, fourth and sixth order, which are related 
to those expansion coefficients, increase with decreasing 
quark mass. 

Alternatively to the quark chemical potentials one can 
introduce chemical potentials for the conserved quantities 
baryon number B, electric charge Q and strangeness S 
(fJ,B,Q,s)i which are related to ix u ,d,s via 



1 2 

Mu — g Ms + 3^Q' 

1 1 

[hi = gM-B - gMQ, 

1 1 



Ms- 



(G) 
(7) 
(8) 



By means of these relations the coefficients cf :?;' s of the 
pressure expansion in terms of hb,q,s are easily obtained, 
in analogy to Eq. [5] 



^ = E«5rm(f)'(fy'(f)*. 



(9) 



For the rest of this article we will restrict ourselves to the 
case of fiQ = [is = 0, thus we will suppress in the following 
the indices that are related to those chemical potentials. 
From the pressure we immediately obtain the baryon num- 
ber density tib, which is given by the derivative of p/T 4 
with respect to the baryon chemical potential /ib and can 
be expressed in term of the expansion coefficients c^, we 
have 

ncf(T)(^) n ~\ (10) 



n B 
T 3 



OO 
71=2 



Using standard thermodynamic relations we can also 
calculate the expansion coefficients of the trace anomaly 
or equivalently the difference between energy density 
and three times the pressure, 



3p 



rp4 



n=0 
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where the expansion coefficients are given by 



c B AT)=T 



dcg(T) 
dT ■ 



(11) 
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Fig. 4. Taylor coefficients of the pressure in term of the up-quark chemical potential. Results are obtained with the p4fat3 
action on iV T = 4 (full) and N T = 6 (open symbols) lattices. We compare preliminary results of (2+l)-flavor a pion mass of 
m-n ~ 220 MeV to previous results of 2-flavor simulations with a corresponding pion mass of m v w 770 |16| . 




Combining Eqs.OHTJ and fT2l we then obtain the Taylor 
expansions for the energy and entropy densities [19) 

oo 

^=E( 3 *)+ c t(T)) {^y 

n=0 

-£«.(£)". M 

n=0 

a _ e +p- \XBriB 

oo 

= £((4-n)c£(T)+i£(T)) (^f) 

n=0 
n=0 

At present, we calculate the expansion coefficients from 
the coefficients , in accordance with Eq. Q21 by perform- 
ing the T derivative numerically, which introduces a small 
systematic error. 



In Fig. [5] we show the second, fourth and sixth order 
expansion coefficients of the pressure, energy density and 
entropy density as given in Eqs. [5] and 1141 obtained with 
the p4fat3 action. Full symbols are from N T = 4 lattices, 
while the open symbols denote results from N T — 6 lat- 
tices. We again find small cut-off effects, however, higher 
order derivatives of pressure, energy density and entropy 
density with respect to \xb are still very preliminary, as the 
error bars are large. This is especially true for the results 
from N T — 6 lattices. Nevertheless, the overall pattern of 
the coefficients is in agreement with expectations based 
on an analysis of the singular behavior of the free energy, 
making use of an appropriate scaling Ansatz. 

We find that the magnitude of the coefficients is de- 
creasing drastically with increasing order, for all analyzed 
temperatures. Thus an approximation of the equation of 
state for small baryon chemical potential by means of a 
fourth or sixth order expansion seems to be justified. In 
general, an analysis of the radius of convergence of such 
a Taylor series is of great interest for an analysis of the 
QCD phase diagram, since the radius of convergence is 
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Fig. 6. On the left panel we show Isentropic trajectories in the (T — /is)-diagram, corresponding to s/ub = 300,45,30 
respectively. Open symbols are form N T — 6, obtained by a 4th order Taylor expansion of the pressure. Filled symbols are from 
N T = 4 calculations. We also show freeze-out data, as well as a parameterization of the freeze-out curve from |21] , Corresponding 
free gas limits are (j,b/T — 0.21, 1.40, 2.13 respectively and are indicated by solid lines. On the right panel we plot the the ratio 
of pressure and energy density along those trajectories. 



bounded by the location of the QCD critical point as well 
as by any first order phase transition line. 



6 The isentropic equation of state 

By using the Taylor expansion coefficients of the baryon 
number (Eq. 1 1 0[) and entropy density (Eq. 1 14j) . we can 
compute the ratio of entropy per baryon number as func- 
tion of T and (1b ■ Solving numerically for a constant ratio 
of entropy per baryon number, s/ns, we determine isen- 
tropic trajectories in the (T, /xs)-plane. These trajectories 
are relevant for the description of matter created in rela- 
tivistic heavy ion collisions. After equilibration the dense 
medium created in such a collision will expand along lines 
of constant entropy per baryon. It then is of interest to 
calculate thermodynamic quantities along such isentropic 
lines. 

We find that isentropic expansion at high temperature 
is well represented by lines of constant [J-b/T down to 
temperatures close to the transition, T ~ 1.2Tq. In the low 
temperature regime we observe a bending of the isentropic 
lines in accordance with the expected asymptotic low tem- 
perature behavior. The isentropic expansion lines for mat- 
ter created at SPS correspond to s/ns — 45 while the 
isentropes at RHIC correspond to s/ns — 300. The en- 
ergy range of the AGS which also corresponds to an energy 
range relevant for future experiments at FAIR/Darmstadt 
is well described by s/ns — 30. These lines are shown 
in Fig. [S] (left) together with data points characterizing 
the chemical freeze-out of hadrons measured at AGS, SPS 
and RHIC energies. These data points have been obtained 
by comparing experimental results for yields of various 
hadron species with hadron abundances in a resonance 
gas 20, 21j. The solid curve shows a phenomenological pa- 
rameterization of these freeze-out data [21] . In general our 



findings for lines of constant s/ns are in good agreement 
with phenomenological model calculations that are based 
on combinations of ideal gas and resonance gas equations 
of state at high and low temperature, respectively [22"ll23|. 

Results shown in Fig. [6] are based on a fourth order ex- 
pansion of the pressure. We find, however, that the trunca- 
tion error is small, i.e. the results change only little when 
we consider also the sixth order term in /is - In accordance 
with the good convergence of our results, we find, that all 
trajectories shown in Fig. El (left) are well within the ra- 
dius of convergence of the Taylor series. At present we 
estimate the radius of convergence of the pressure series 
to (/! B /T) crit >2.7. The cut-off effects can be estimated by 
comparing open and full symbols. 

We now proceed and calculate energy density and pres- 
sure on lines of constant entropy per baryon number using 
our Taylor expansion results up to 0(/j, b ). We find that 
both quantities obtain corrections of about 10% at AGS 
(FAIR) energies (s/ns = 30) and high temperatures. The 
dependence of e and p on s/ub cancels to a large extent 
in the ratio p/e, which is most relevant for the analysis of 
the hydrodynamic expansion of dense matter. This may 
be seen by considering the leading 0(fj, B ) correction, 



1 eg - 3p 
3 eo 



C-2 



e - 3p 



£2 



ffi 



(15) 

In Fig. El (right) we show p/e as function energy density 
along our three isentropic trajectories. The softest point of 
the equation of state is found to be (p/e) m in — 0.07 — 0.09, 
for N T = 4 and 6 respectively. Within our current numer- 
ical accuracy it is independent of s/ns ■ Similar results for 
the asqtad action have been obtained in [18] , However, as 
our data is preliminary, the analysis clearly suffers from 
poor statistics, which is in particular true for our iV r = 6 
results. 
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7 Conclusions 

We have presented results on the equation of state on 
lattices of N T = 4,6 [T|IT4] and 8 with two differ- 
ent kinds of improved staggered fermions. Our masses 
have been kept constant in physical units and are cho- 
sen such that we have a physical strange quark mass (m s ) 
and 2 light quarks with a mass of mi = 0.1m s . We also 
presented some preliminary results with physical quark 
masses mi = 0.05m s . We find that our two actions lead to 
a consistent picture of the thermodynamics of QCD and 
find in particular for the N T = 8 results only small cut-off 
effects. We have calculated the equation of state as well 
as the velocity of sound and find the softest point of the 
equation of state to be (p/e) m in — 0.09 at energy densities 
of 1 - 2 GeV/fm 3 . 

Furthermore, we calculated corrections to the equation 
of state arising from a non-zero baryon chemical potential, 
by means of a Taylor expansion of the pressure. Within 
this framework we calculated the isentropic equation of 
state along lines of constant entropy per baryon number 
{s/n B ) for RHIC, SPS and AGS (FAIR) energies. Within 
our current, preliminary, analysis we find the softest point 
of the equation of state to be independent of s/ub- 
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